Emergence of Plasmaronic Structure in the Near Field Optical Response of Graphene 
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The finite momentum optical response cr{q, u) of grapfiene can be probed with the innovative 
technique of infrared nanoscopy where mid-infrared radiation is confined by an atomic force micro- 
scope cantilever tip. In contrast to conventional g ~ optical absorption which primarily involves 
Dirac fermions with momentum near the Fermi momentum, k ~ kp, for finite q, G{q,ui) has the 
potential to provide information on many body renormalizations and collective phenomena which 
have been found at small k < kp near the Dirac point in electron-doped graphene. For electron- 
electron interactions, the low energy excitation spectrum characterizing the incoherent part of the 
quasiparticle spectral function of Dirac electrons with k r^ kp consists of a flat, small amplitude 
background which scales with chemical potential and Fermi momentum. However, probing of the 
states with k near fc = will reveal plasmarons, a collective state of a charge carrier and a plas- 
mon. These collective modes in graphene have recently been seen in angle-resolved photoemission 
spectroscopy and here we describe how they manifest in near field optics. 

PACS numbers: 78.67. Wj,73.22.Pr,73.20.Mf 



INTRODUCTION 

Since its discovery, graphene continues to provide a 
rich trove of exotic phenomena arising out of novel charge 
dynamics vifhich map onto a Hamiltonian for massless 
Dirac fermions, providing a low energy bandstructure 
of two linear bands which touch and cross at a so- 
called Dirac point. Using angular resolved photoemis- 
sion spectroscopy (ARPES) Bostwick at al. [1] have re- 
cently observed, in samples of graphene grown epitax- 
ially on a H-SiC substrate, that the Dirac point splits 
into two at momentum fc = 0. A plasmaron ring is ob- 
served between the two Dirac points which has its ori- 
gin in electron-electron interactions (EEI) [2-1^1 and con- 
sequently is not part of the bare band description. In 
metals, the long range nature of the Coulomb interac- 
tions provides the mechanism for the formation of long- 
lived collective charge density oscillations or plasmons. 
These also exist in graphene [5|, |6j although with an al- 
tered dispersion relation, i.e., in the long wavelength 
limit {q — >■ 0) there is a square root dependence on q. As 
recognized by Lundquist [7| , charged quasiparticles can 
interact with the plasmons to form a composite state, 
termed a plasmaron. It is these composite excitations 
that provide an understanding of the ARPES data in the 
region around A; ~ near the split Dirac point. Other 
evidence of plasmarons in graphene might be found in 
the tunneling conductance. The observedfsl, |9| lifting 
of the node in the bare band density of quasiparticle 
states provides evidence of the presence of many-body 
effects 1Ch13| , moreover plasmaronic structures in this 
region, which scale with chemical potential, 1^ |13| , may 
underlie the EEI-type features identified [9[ in tunneling 
experiments. Aside from graphene, the elusive plasmaron 



has been observed in other systems [lj| such as elemental 
Bi semimetals. la ll6| 



Optical spectroscopy provides useful information on 
charge carrier dynamics and has widely been used to 
probe the nature of the Dirac quasiparticles in graphene. 
Many experimental results are reviewed in the paper 
of Orlita and Potemski.'l7| While most of the observa- 
tions can be understood qualitatively and even quanti- 
tatively with bare bands, some cannot. A prominent 
example is the absorption spectrum of graphene as de- 
scribed by the real part of its AC conductivity at 9 = 0, 
(j{q = OjCj). Theory predicts a Drude peak centered 
about zero photon energy (w) followed by a region of 
near zero absorption before a rapid rise to the universal 
value of cTn = 7re^/2/i at w = 2/x, where ^ is the chemical 
potential. [18l42d | This background remains nearly con- 
stant over a large energy range above 2/x. The Drude 
piece comes from intraband absorption [see lefthand side 
of Fig. [Tfb)] and its width is set by the residual scatter- 
ing rate 7 assumed to be much smaller than 2/i. The 
interband optical transitions, which provide the univer- 
sal background, begin only at uj — 2 fj, because of Pauli 
blocking, i.e., Pauli exclusion blocks vertical transitions 
for photon absorption below this energy due to an oc- 
cupied final state for k < kp [Fig. [ijb), q = 0]. One 
aspect of this bare band picture is not confirmed in exper- 
iment. The conductivity never falls below about (To/3 in 
the Pauli-blocked region of the spectrum where it would 
be expected to be almost zero.[21| This has been widely 
interpreted as evidence for the importance of Coulomb 
correlations and the electron-phonon interaction which 
provides phonon-assisted Holstein sidebands ll|, |22| to 
the main Drude peak. A second contribution comes from 
equivalent sidebands associated with the rise of the in- 
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FIG. 1: (Color online) (a) Schematic of a graphene sheet on a 
substrate with an AFM tip used to confine the incident radi- 
ation to a nanoscale region around the tip as in Fei et al.[23]. 
(b) Schematic of intraband and interband transitions in the 
Dirac cone bandstructure of graphene. Transitions from occu- 
pied (shaded) to unoccupied states for g = photon absorp- 
tion are shown on the left with the heavy shading indicating 
the Pauli-blocked region. For a finite q transfer, new optical 
absorption processes (shown at right) are possible. For ex- 
ample, a particle at k = kp — q may be excited to k = kp 
with finite q transfer. The Pauli-blocked region is reduced as 
q increases and the plasmaron region around fc = can be 
unveiled and probed by optics. 



terband optical transitions at w = 2/i.|ll| 

In this letter, we address the issue of how plasmaronic 
structures present themselves in the optical conductivity 
of graphene. Conventional experiments measure the q ^ 
0, long wavelength limit, of the AC conductivity a{q,uj). 
For this case the intraband transitions mainly involve 
the quasiparticle states at and near momentum k = kp 
in the top Dirac cone, while the interband transitions 
involve k > kp for both upper and lower Dirac cones 
[lefthand side of Fig.[ljb)]. Consequently, the plasmaron 
region found around k ~ (as we have described) is 
not directly probed and one does not expect to see a 
significant signature of the plasmarons in the g ~ AC 
conductivity. However, this is not the case for the finite 
q conductivity. 

There have been significant advances in infrared 
nanoscopy which effectively allow for the measurement of 
the finite q optical conductivity cr(q,a;) as a function of 
(jj. A cantilever tip of an atomic force microscope (AFM) 
operating in a tapping mode (as depicted in Fig. [T{a)) 
allows for the acquisition of near field data containing 
information on the finite q conductivity. In the paper by 
Fei et al. 



23| a confinement region of « 30 nm is achieved 



super sharp silicon tip might reduce the confinement to 
10 nm and tips based on carbon nano-tubes to perhaps 
1 nm. As an example, one sees from the ARPES data 
presented in Fig. 2(b) of Bostwick et al.[l| [and shown 
theoretically here in Fig. [H^a)] that the Fermi momen- 
tum, kp ^ 60 X fO^'^A^^, for ^ of roughly 600 meV. 
The plasmaron region falls between 500 and 800 mcV 
and corresponds to a momentum region about fc = of 
roughly 20 x lO^'^A^^. To access this region a g of or- 
der 40 X lO^'^A^^ is required which is beyond the reach 
of the tip used by Fei et al.[23] but in the range that 
could be achieved with a carbon nanotube tip. How- 
ever, reducing the doping to obtain a chemical potential 
of 200 meV would bring the presently available q close 
to the plasmaron region. Indeed, previous g = opti- 
cal experiments [2 15 were performed on samples on Si02 
substrates with n ^ 150 — 300 meV. The scaling of the 
plasmaron feature with /i implies that the current range 
of q values might already be able to access the plasmaron 
ring region. 

The Kubo formula for the real part of the finite q op- 
tical conductivity is given in terms of the universal back- 
ground conductivity ctq by 
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where A^{k,uj) is the spectral function of the Dirac 
fermions and Fss'{(j)) — ^[1 + ss' cos (j)]. Here (j) — 9 + (3, 
where /? = arctan( ^^'"osg ) and 9 is the angle between 
fc and q. Note here that (j> differs from that associ- 
ated with the polarization function, H(q,w), wherein 
(f) — 9 — p. While for bare bands with q — >■ 0, a{q,uj) = 
Je^a;H(q,cj)/q^ as seen in Ref. [2J], at finite q and with 
interactions this relation does not hold. For graphene, 
the spectral function can be written in terms of the self- 
energy for the s = ± cones, Ss(fc,w), as 
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providing a distribution of q values which peaks around 
q = 3.4 X 10^ cm~^. As described by these authors, a 



(2) 

where the bare band energies are given by e| — 
svpk — fi. In this work we use the standard GqW- 
RPA approximation [12| for the self-energy arising from 
Coulomb interactions (EEI) . We have set H = 1 and will 
also set the Fermi velocity vp = 1 when discussing en- 
ergy transfer. In Fig. [TJb) we show the Dirac cones for 
finite chemical potential fi. The plasmaron region which 
we would like to access falls near fc = and extends to 
fc = fc"*" as shown in Fig.[5Ja). A crossing occurs between 
a plasmaron band and a Dirac quasiparticle band in this 
region. The plasmaron structure in Fig.[2]scales with the 
value of the chemical potential /x, and Fermi momentum, 
kp, and applies to any /x. Returning to Fig. [T][b) (right), 
finite q optical transitions are shown in which a hole, de- 
picted as an open circle is created by the light with the 



excited final state electron at k = kp, shown as a solid 
(red) circle. It is clear from this picture that these new 
intraband and interband transitions sample states in the 
Dirac cone down to k ^ kp — q which approaches the 
Dirac point at fc = as g approaches the Fermi momen- 
tum. To be close to, or even inside of, the plasmaron ring 
region, a momentum oi q ^ kp — k^ is needed. The over- 
all effect of finite q probing is to shrink the Pauli-blocked 
region towards the Dirac point, thus uncovering the re- 
gion of fc-space where plasmarons are most prominent. 

Plasmaronic modes and dressed Dirac quasiparticles 
are found in solutions of the equation 



Rel]^(fc,w) -e| ==0. 



(3) 



The former are introduced through the interactions and 
represent new collective modes (or scattering resonance) 
of a charge carrier plus a plasmon. The contours along 
which Eq. ([3]) is satisfied are indicated by the compli- 
cated structure shown in Fig. [2Kb) with the upper conical 
region representing the renormalized Dirac quasiparticle 
branch and the lower cone, the plasmaronic one.pj This 
has been probed in detail in ARPES work of Bostwick et 
al.Ql who also confirm that GqW-RPA theory provides 
a good description of the spectral function of the Dirac 
fermions in graphene. A color plot of our own results for 
the spectral density, A(fc,w), is shown in Fig. [2Ka) with 
emphasis on the region near k — Q. One of the moti- 
vations of this work is to encourage confirmation of the 
ARPES results using optical methods. 

Another representation of the plasmaron structure is 
shown in the bottom frame of Fig. [2] where we provide 
in the inset a three dimensional plot of the Dirac carrier 
spectral density. Around k = kp we see a sharp quasipar- 
ticle peak with a very diffuse incoherent sideband extend- 
ing to large energy but always flat and small in amplitude 
for w > — 2/z, as is made clearer in the solid (black) curve 
of Fig. EKc). Below ui ~ —2^ there is a plasmaron peak 
around — 2.8/i with a small sideband from Dirac quasi- 
particles but this appears at large negative energies well 
within the region of the universal background in optics 
therefore not of primary interest here. As k is reduced the 
plasmaron structure becomes more prominent in the en- 
ergy window of interest and the dash-dotted (blue) curve 
for k = k^ shows three peaks: the main quasiparticle 
peak near — 0.9/i, a plasmaron peak at ^ — 1.8/i and a 
peak related to the plasmaron ring (the crossing of the 
plasmaron and main quasiparticle bands) in between. By 
/c = 0, two peaks of almost equal magnitude remain in the 
dashed (red) curve, the upper one marks the quasiparticle 
Dirac point crossing and the lower, the plasmaron Dirac 
point. It is this rich structure from ^ k^ towards fc = 
which would be exciting to access with optics. However, 
there does not exist a sharp transition in the evolution of 
A{k,uj) as a function of fc at fc = fc+, but rather the tran- 
sition is gradual and important structure representative 
of electron-electron correlations comes into play before 
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FIG. 2: (Color online) (a) Color map of the Dirac fermion 
spectral function A{k,uj) including EEL (b) Solutions of 
Eq. (|3j for quasiparticle energies dressed by EEI in GoW-RPA 
theory. A complex set of additional solutions to Eq. ((3| ap- 
pear between the split Dirac point and extend into the lower 
cone, which correspond to new collective states called plas- 
marons. (c) A{k, uj) vs u) for k = kp, fc"*", and 0. Inset: three 
dimensional view of A{k, ui) showing constant momentum cuts 
(solid color lines) for these three k values. 



k — k^ and at energies larger than ~ — 2.0/i and this is 
also of great interest. 

In Fig. [31 we present our results for the finite q optical 
conductivity a{q,u]) vs w for five values of q/kp = 0.0, 
0.1, 0.4, 0.8 and 1.0. The curves are for the case of no 
EEI. In the numerical work, we have used a constant 
residual scattering rate, 7 of 0.015. The solid black curve 
for gr = is well known. It has a Drude peak centered 
about u — Q which decays rapidly according to the value 
of 7. The region of near zero conductivity is followed with 
a sharp rise towards the universal value ctq at cj > 2/i. As 
q is increased a{q, u) shows a new peak at a; = g (recall, 
vp — \ here) corresponding to an average of displaced 
oscillators with a; = ^k+a ^ ^k which has a maximum 
value of q. This is followed by a drop in the conduc- 
tivity which becomes progressively less as q is increased. 
There is a point of inflection at uj — 2^ — q, the energy 
at which the interband transitions begin [see Fig. [IJc)]. 
Finite q provides absorption in the region which would be 
Pauli-blockcd for g = and this increases with increas- 
ing q. Another important energy scale is 2/i -I- g above 
which the universal plateau shows no depletion and be- 
low which there is some reduction. All the curves meet at 
a single point w = 2/x and this crossing could be used to 
identify the chemical potential. The first peaks in each 
of the curves for a{q, uj) vs w is correlated directly with 
the quasiparticle peak in the spectral density, A(k — q, uj) 
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FIG. 3: (Color online) (a) The real part of the finite mo- 
mentum optical conductivity (j{q, uj) as a function of uj/fj, for 
q/kp — 0, 0.1, 0.4, 0.8 and 1.0 for bare bands. Prominent 
quasiparticle peaks at u — q are identified. The interband 
edge which rises sharply at oj = 2/i for q — becomes pro- 
gressively less prominent with increasing q and shifts to 2fi — q. 
Another feature is seen at u — 2fi + q below which the univer- 
sal background is depleted for finite q. (b) The finite q optical 
conductivity a{q,oj) for g=0.4, compared with the spectral 
function A(k, u) for k — kp — q and k — kp + q. 



(dashed blue curve) of Fig.[3ljb). The spectral density is 
the sum of contributions from the top and bottom Dirac 
cones so that there are two peaks in A{k,uj). The sec- 
ond peak at high frequency correlates with the onset of 
the interband transitions. The dashed-dotted curve is 
A{k + q, to) and this peak is at w = 2/x -|- g where there is 
a kink in the conductivity. These curves are to be com- 
pared with those in Fig.|3l^a) which have the same format 
as the bare band results but include the EEI in Eq. ^. 
Each curve displays a clear quasiparticle peak displaced 
from the bare case as a result of the appearance of the 
real part of the self-energy in Eq. ([3]) . More importantly, 
new structure (identified by arrows) appears in a{q,uj), 
not seen in the bare band case and are due to correla- 
tion. While such effects are, in principle, also present in 
the q = case, they provide only a small modulations [25| 
on the shape of the intraband Drude peak which depends 
only on electrons with k = kp- As seen in the inset of 
Fig. [2Ic), correlation sidebands are very weak for these 
electrons. But EEI structure is sufficiently sharp to show 
clearly in the q/kp = 0.4 curve. This comes from the 
plasmaron sideband shown in Fig. [2ja) which grows in 
amplitude and sharpens as k is reduced from k = kp 
to k = 0. The optical conductivity for q/kp = 1 in- 
deed samples the bands converging to the two split Dirac 
points at /c = and in this case the first peak around 
w ~ 1.1/x is the Dirac point from the main quasiparticle 



FIG. 4: (Color online) (a) The real part of a{q,uj) for 
q/kp^O-O, 0.1, 0.4, 0.8, and 1.0 including EEI. Quasiparti- 
cle peaks are seen in each curve but displaced from their bare 
band positions and these are followed by new plasmaronic 
shoulders indicated by arrows. Interactions increase the fill- 
ing of the Pauli-blocked region as emphasized in frame (b) 
where interacting and non-interacting cases are compared. 



bands (shifted slightly from uj = fi hy the renormaliza- 
tions) and the second hump around w ~ 1.6/i is due to 
the plasmaron band Dirac crossing. To make our point 
clear we plot the positions of the peaks in Fig. HJ^a) for 
each q on the band structure of Fig. [^fa) using w — ?> — w 
and |fc| = \kp — q\. These are shown as open triangles 
(circles) and track the plasmaron (quasiparticle) band. 
This demonstrates how measurements from optics could 
be used to trace the renormalized band structure. In 
Fig. mjb), the dashed green curve is reproduced from 
Fig.[3lja) and is shown for easy comparison with the solid 
green curve for the same q/kp = 0.8 which includes EEI. 
Not only is the quasiparticle peak shifted and a plas- 
maron shoulder seen, but interactions have filled in the 
large dip present in the case of the bare bands 

In summary, for the bare bands the intraband optical 
transitions, which provided a Drude peak in a{q,(jj) for 
q = 0, are now, for finite q values, distributed over an en- 
ergy range with the dominant peak moved to u — q. This 
fills the Pauli blocked region of the g = case. Further- 
more, the interband transitions now start at uj = 2^, — q 
rather than 2//, and the universal background ctq shows 
some depletion up to 2/x -I- q. Electron-electron interac- 
tions shift the bare band peak atuj — qio higher energies, 
refiecting quasiparticle energy renormalizations. They 
also provide further filling of the Pauli-blocked region 
and show easily identifiable new shoulder-like structures 
which have their origin in the plasmaronic side bands of 
the Dirac fermion spectral density, A{k,uj). The signa- 
ture of these modes becomes progressively more promi- 



nent as the region about fc = is sampled. 
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